Non-trivial fixed point structure of the two-dimensional ±J 3-state Potts 

ferromagnet/spin glass 
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The fixed point structure of the 2D 3-state random-bond Potts model with a bimodal (±J) dis- 
tribution of couplings is for the first time fully determined using numerical renormalization group 
techniques. Apart from the pure and T = critical fixed points, two other non-trivial fixed points 
are found. One is the critical fixed point for the random-bond, but unfrustrated, ferromagnet. The 
other is a bicritical fixed point analogous to the bicritical Nishimori fixed point found in the random- 
bond frustrated Ising model. Estimates of the associated critical exponents are given for the various 
fixed points of the random-bond Potts model. 



Over the past decades the study of random systems 
has attracted considerable attention. In particular the 
possibility of a spin-glass phase has been studied in great 
detail, as a result of which it now seems well established 
that no finite temperature equilibrium phase transition 
to a spin-glass phase occurs in two dimensions [1]. How- 
ever, even in the absence of a glassy phase a rich fixed 
point structure can occur. In particular, the transition to 
an ordered phase in the presence of disorder can be con- 
trolled by new random fixed points. Recently, it has been 
suggested that while pure systems can display critical be- 
havior in many different universality classes, the number 
of universality classes may be more limited when disor- 
der is present [2]. In addition, current progress in the un- 
derstanding of quantum phase transitions in disordered 
systems as well as other complex systems has shown that 
many of these models are in the same universality class as 
certain disordered classical statistical mechanical mod- 
els [3]. It is therefore highly desirable to undertake a 
careful investigation of the complete fixed point structure 
of some key models, determining all stable and unstable 
fixed points. While considerable progress towards the 
understanding of random fixed points has been made us- 
ing analytical techniques [4-7] and series expansions [8] , 
it would seem that a first step towards such a classifi- 
cation of the fixed point structure would have to come 
from numerical work that does not rely on perturbative 
methods and where frustration can be taken into account 
satisfactorily. McMillan [9] investigated the fixed point 
structure of a two-dimensional, frustrated random-bond 
Ising model using the domain wall renormalization group 
(DWRG) method. In this case a bicritical point occurs, 
separating the ferromagnetic-to-disordered critical phase 
boundary into two parts. The part at higher tempera- 
ture and weaker disorder is governed by the pure criti- 
cal fixed point, while the lower-temperature portion of 
the phase boundary is governed by the zero-temperature 
fixed point that also governs the ferromagnetic-to-spin 
glass transition at zero temperature. For the random 



Ising model this bicritical point occurs on the so-called 
Nishimori line [10], where there is a special gauge sym- 
metry, and the ferromagnetic and spin-glass correlations 
and susceptibilities are identical [8]. 

In the present paper we apply a variant of the renor- 
malization group (RG) method of McMillan [9] to the 
3-state Potts model and determine all the fixed points. 
For this model it is known from the Harris criterion [11] 
that the pure fixed point is unstable against disorder and 
the simplest possible scenario for the RG flow would be a 
flow out from the pure fixed point that goes to the zero- 
temperature fixed point separating the ferromagnetic and 
spin-glass phases. However, we show that this is not the 
case and that two additional non-trivial random fixed 
points occur at nonzero disorder and nonzero tempera- 
ture. In agreement with analytical work using only fer- 
romagnetic (unfrustrated) disorder [4,6], an unfrustrated 
random critical fixed point is found. However, in addi- 
tion to this critical fixed point a bicritical fixed point 
also occurs in this model, even though the model does 
not have the gauge symmetry that is used to define the 
Nishimori line in the Ising model. We believe that this 
fixed point structure is of a rather general nature, in the 
sense that it should be similar for many other models for 
which disorder is relevant. 

The Hamiltonian we use is: 



H = - ) Jijh(ni,n,j), n,- = 1 . 



.q, 



(1) 



<ho> 



with h(ni, rij) — cos(2?r(nj — rij)/q), and q — 2 (Ising) or 
q = 3 (Potts). The sum is over all nearest-neighbor pairs 
of sites on a square lattice. We consider the situation 
where the bonds Jy in Eq. 1 are distributed randomly, 
and given by a quenched biased bimodal probability dis- 
tribution: V{J) = x8(J - 1) + (1 - x)S(J + 1). As an 
application of this type of model, it has been suggested 
that the orientational freezing in molecular glasses, such 
as N2-Ar and KBr-KCN, can be partially described by a 
three-dimensional 3-state Potts spin glass [12]. 
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The way we apply the DWRG is slightly novel and we 
therefore begin by reviewing the method and the scal- 
ing ideas behind our approach. The DWRG as pro- 
posed by McMillan [9] estimates the singular part of 
free energy by calculating the domain wall energy [AF], 
and its standard deviation, cr(AF), for systems of size 
L x L. Here [•] denotes disorder averaging. We calcu- 
late the total free energy difference between periodic (p) 
and twisted (t) boundary conditions, SF m = F t — F p , 
for a long cylinder of size L x M where M = Lm, and 
m is a large integer running up to of order 10 5 — 10 6 . 
The proper generalization of twisted boundary conditions 
to q > 2 is to change the local Hamiltonian across the 
boundary in the L-direction as follows [13]: h(n\,ni) — 
cos(27r((ni -nt + l)/</)- This change is done at a seam 
going along the entire length of the cylinder. AF for each 
LxL subsystem is simply defined as AF = <5F m — <5F m _i. 
The free energies are evaluated exactly using standard 
transfer matrix techniques [14] and the only errors in the 
calculation therefore stem from the incomplete disorder 
averaging. 

Pure systems. In the absence of disorder it is known 
from hyperscaling that the singular part of the free en- 
ergy density scales as the inverse of a correlation volume 



fs _C 

hT c e' 



(2) 



Since f s ~ AF/L d it follows that an appropriate finite- 
size scaling form for AF is 



AF 



(3) 



with g a universal function and 5 = \T — T c \. Hence, the 
critical point, T c , can be located by standard methods, 
i.e. by tracing AF as a function of T for several different 
system sizes and locating the point where the lines cross. 
Here A is the universal amplitude for the spin stiffness of 
the system. The universality of A has been investigated 
extensively at finite temperature transitions [13,15-18], 
and is known exactly for the pure q — 2,3,4 Potts models 
as a function of the aspect ratio s = 1/m [13]. In the 
limit s — > it follows from conformal invariance that 
A = ttt] [16], where r\ is the magnetic critical exponent. 
Since we use an aspect ratio of s = 1/m, essentially zero, 
we see that = nn. Hence, AF directly measures a 
bulk critical exponent. 

Disordered systems. Just as AF/ksT is a universal 
amplitude at the pure fixed point it is natural to assume 
that any random fixed point occurring at finite, nonzero 
temperature and disorder will be characterized by a uni- 
versal distribution of AF/fceT. In particular, the mean 
[AFj/fceT and the standard deviation a(AF)/k B T are 
then universal [19,20]. The generalization of the finite 
size scaling relation Eq. (3) to include disorder is there- 
fore: 



[AF] 
hT 



(4) 



valid close to any finite temperature fixed point X 
at (x*,T*), with gx a universal scaling function that 
is, however, different for each distinct fixed point X. 
cr(AF) /ksT obeys a similar finite-size scaling form. Here 
6 is a linear combination of Air = \x — x* |, AT = \T— T*\ , 
that is the eigenvector of the RG flow that moves one 
away from the phase boundary; 1/v is the corresponding 
leading eigenvalue, e is another linear combination that 
is the next subrelevant (A2 > 0) or irrelevant (A2 < 0) 
eigendirection of the RG flow that is tangent to the phase 
boundary. Thus we can estimate the fixed points (x*,T*) 
as the points in the phase diagram where ([AF], cr(AF)) 
are independent of L. We obtain a picture of the RG flow 
by measuring [AF] and cr(AF) at various (x, T) points 
and seeing how they vary as L is increased at constant x 
and T. 

Zero temperature. If the transition is controlled by 
a zero temperature fixed point Z, then hyperscaling, 
Eq. (2), is not valid, and we cannot use Eq. (3) as our 
starting point. The appropriate finite-size scaling form is 
instead 



[AF}=L e g z (6L 1 /»,eL- e ), 



(5) 



with 9 > 0. Hence, [AF] is not universal along the part 
of the phase boundary where the flow is controlled by a 
zero-temperature fixed point. Instead, [AF] grows with 
increasing L as L e at the critical point. At small tem- 
peratures we expect that S ~ \x — x c \ and T ~ e, i.e. 
that the disorder x is the relevant flow direction. A pos- 
itive 9 implies that the flow is into the zero-temperature 
fixed point and hence that temperature is an irrelevant 
variable. 

We now proceed as follows: In order to accommo- 
date both finite temperature and zero temperature fixed 
points in the same plot we consider the two quantities 
r = cr(AF)/[AF] and / = k b T/[AF]. Any fixed point 
should be characterized by universal values of these quan- 
tities (r*,/*). In particular we know that for the pure 
fixed point we have (r*,f*) = (0, (ttt])^ 1 ). A ferro- 
magnetic phase will be characterized by (r, /) — ► (0, 0) 
whereas in the paramagnetic phase [AF] — ► 0. we can 
now directly determine the complete fixed point struc- 
ture by calculating (r(L),/(L)) for different values of L 
on a grid of (x, T) close to the phase boundary. By con- 
necting points, (r(L),/(L)) corresponding to L and L' 
calculated at the same (x,T), with an arrow terminat- 
ing at the larger of L, L', the flow becomes clearly visible 
and "fixed points" (r*, /*) appear where both / and r are 
independent of L. The linear flow around a fixed point 
should be characterized by 



f(L') 
r(L') 



f* 
r* 



r 

r* 



(6) 
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Through a least-square fit to this equation for a num- 
ber of points close to the fixed point the matrix M as 
well as (r*,/*) can be determined and consequently the 
eigenvalues, Ai = A2, and eigenvectors of the linear 
flow can be determined. The magnetic exponent 77 is es- 
timated by exact calculations of [< m 2 >] ~ L 2 ~ r ', the 
square of the magnetization per site, directly at the fixed 
point. 
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figure 1: The flow close to the pure fixed point for 
the 2D Ising model with bimodal ferromagnetic disor- 
der V(Jij) = (1 - x)d(Jij - 1) + x5(Jij - 0.3). P is the 
exactly known pure fixed point. In most cases 1 million 
L x L blocks were measured. 

In order to establish that our approach is sensitive even 
to marginal flow we first consider the two-dimensional 
Ising model with ferromagnetic disorder. In this case it 
seems well established [6,4,5] that this type of disorder 
is marginally irrelevant at the pure critical point. Hence 
the flow should be towards the pure critical fixed point P. 
Our results for this case are shown in Fig. 1 for a bimodal 
distribution of couplings of strength J and 0.3 J. The 
arrows connect results for L = 6 and L' = 8. It is clear 
that the flow in this case indeed is towards P, located at 
(0, 4/7r), in agreement with the theoretical studies [6,4,5]. 

We now turn to a discussion of our main results for the 
3-state Potts model. The flow diagram is shown in Fig. 2 
with L = 6 and V = 8. Rather clearly, four different 
fixed points are visible along the phase boundary of the 
ferromagnetic phase. They are: the pure critical fixed 
point P, the random critical fixed point S, the bicriti- 
cal fixed point N and the zero-temperature fixed point 
Z. Close to the latter fixed point our calculations have 
been performed exactly at zero temperature. The RG 
flow along the phase boundary is away from P and N 
and towards S and Z. Also indicated (bold arrows) in 
Fig. 2 are the numerically determined RG eigenvectors at 
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figure 2: The flow diagram for the q=3 Potts model. P 
is the exactly known pure critical fixed point, S the ran- 
dom critical fixed point, N the analogue of the Nishimori 
bicritical fixed point, and S is the zero temperature fixed 
point. The bold arrows indicate the numerically deter- 
mined RG eigenvectors. The PM and FM indicate the 
paramagnetic and ferromagnetic phase. In most cases 
200,000 L x L blocks were measured. 



the fixed points S and N. At the stable fixed point (S), 
with approximately (x*,T*) - (0.93,1.16) (r*,/*) = 

(0.68(9), 0.89(3)) we find Ai = 0.94 - 1.04, A 2 .2, 

and hence v = 1/Ai ~ 0.96 — 1.06. Using [< m 2 >] we 
find r\ — 0.23 — 0.28. Both of these exponents are in good 
agreement with analytical results [4,6] as well as Monte 
Carlo work [21], for the q = 3 Potts model with ferromag- 
netic disorder, so S appears to be the same random fixed 
point as in the non-frustrated random ferromagnet q=3 
Potts model. The bicritical fixed point (N) at (x*,T*) ~ 
(0.88,0.68), (r*,/*) = (1.01(3), 0.43(3)) is characterized 
by the eigenvalues Ai = 0.55 - 0.65, A 2 = 0.30 - 0.40 
using the lattice sizes L = 6 and 8. However, at this 
fixed point a standard scaling plot of / approximately 
along the eigendirection of Ai and including L = 10 yields 
1/u ~ 0.75 — 0.8. Hence, we estimate v = 1.28 — 1.36, and 
using [< to 2 >] we find r] = 0.17 - 0.22. Both of these 
exponents are remarkable similar to the ones determined 
at the bicritical point in the 2D random Ising model [8] 
of v = 1.32(8), r] = 0.20. The zero temperature fixed point 
(Z) located at x* ~ 0.88(1), r* ~ 0.96 and we determine 
the exponents v = 1.45 — 1.55, n = 0.18(1), where we for 
the Ising model find x* = 0.89(1), v = 1.35 - 1.45, 77 = 
0.18(1), in agreement with McMillan [9]. At this fixed 
point we also estimate 9 ~ 0.3. Since this exponent is 
positive, temperature is irrelevant at this fixed point. 
Several conclusions can now be drawn for the first 
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time for the randomly frustrated 3-state Potts model 
from the results in Fig. 2. The random critical fixed 
point S has exponents that to within our precision ap- 
pear close to those of the pure Ising model, and in good 
agreement with theoretical predictions [4,6]. However, 
we find cr(AF) ^ at S, which certainly distinguishes 
it from the pure Ising critical fixed point. At S we have 
also determined the fraction of negative AF for the fixed 
point distribution. This fraction, which is a measure of 
the "renormalized" frustration level, is decreasing with 
increasing L at S indicating that the fixed point is the 
same as for unfrustrated ferromagnetic disorder. A new 
bicritical fixed point N appears even though the model 
is not gauge-symmetric [10]. This fixed point has expo- 
nents similar to the bicritical fixed point for the Ising 
model [8]. At the zero-temperature fixed point Z tem- 
perature is irrelevant as it is for the Ising model. The 
critical concentration x c is at this point identical for the 
Ising and q — 3 Potts model to within the numerical 
resolution. Recently, Nishimori [22] has rigorously estab- 
lished that X n°"- Isin s > Rising for any bimodal vector-spin 
model that is invariant under the operation Si — > — S[. 
However, the q = 3 Potts model is not invariant under 
this operation and it is therefore remarkable that we for 
the q = 3 model find that the inequality apparently is 
satisfied as an equality. Finally, it would appear that 
the q = 3 Potts model is slightly reentrant in Fig. 2, as 
occurs for the same model on hierarchical lattices [23]. 
However, here, we believe that this is due to the fact 
that corrections to scaling are more pronounced at low 
temperatures, in particular for a bimodal distribution. It 
seems also possible that the phase diagram could show 
reentrance in (r, /) space but not in (x, T) space. For the 
bimodal Ising model it is known that the phase boundary 
is vertical [24] down from the bicritical point (N) where 
the Nishimori line [10] intercepts the phase boundary. 
Remarkably, this is also what we find in the (x, T) space 
for the q = 3 Potts model. Summarizing, our results indi- 
cate that at N and Z the q = 2 (Ising) and q = 3 random 
bimodal Potts models are remarkably similar. However, 
for ferromagnetic disorder the Ising model does not have 
a stable random fixed point at finite temperature [4,6]. 
The stable random fixed point (S) of the q = 3 model 
has critical exponents very close to the pure q = 2 Ising 
model. In light of these observations we believe that the 
fixed point structure shown in Fig. 2 is rather general 
for two-dimensional disordered models when disorder is 
relevant. 
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